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ABSTRACT 

We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously 
selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey 
North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of > 3 at 100 /im down to 
3.0 mJy, and/or at 160 pun down to 5.7 mJy. The fraction of SMGs seen at 160 /im is higher than 
that at 100 /xm. About 50% of radio-identified SMGs are associated with PACS sources. We find 
a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could 
be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at 
higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMC sample) remain 
unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength 
may imply higher redshifts for them than for the average SMC population (e.g., z > 3 — 4). The 
Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far- 
infrared spectral energy distributions have significantly different shapes than template libraries of 
local infrared galaxies. 

Subject headings: Galaxies: high-redshift — Galaxies: individual (GN13 alias HDF850.4) — Galax- 
ies: starburst — Cosmology: observations — Infrared: Galaxies — Submillimeter: 
galaxies 
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1. INTRODUCTION 

Several hundred dust-enshrouded high-z sources have 
been selected through submm/mm imaging with 
bolometer cameras like SCUBA , LAB OCA, AzTEC, 



MAM BO fe.g.. iBertoldi et al.l 



iDannerbauer et all 200: 



Coppin et al 



Hughes et al 



19981; iPerera et al.' '2008t (Pope et al.l 120061 : ISmail eTa! 
2002; Wcifi et al. 200^. The large beam size in 
the (sub) millimeter (e.g., MAMBO: 11"; SCUBA: 
15"; LABOCA: 19") hampers the identification of 
these so-called Submillimeter Galaxies (SMGs; see 
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for a review iBlain et al.l l2002f ) based on bolometer 
data only. The most obvious choice for obtain- 
ing the subarcsecond accurate positions of the dust 
continuum is (sub)millimeter interferometric contin- 
uum observations ( e.g., Dannerbauer et al. 2002, 200^ 
I Younger et al.ll2007( ). However the slow mapping speed 
of (sub)millimeter interferometers does not allow us to 
study a large number of sources. The most suitable tool 
for counterpart identification is interferometric observa- 
tions at radio wavelengths. About 50% — 80% of these 
submm/mm sources have been identified, mainly based 
on radio observations, and the peak of the redshift dis- 
tribution of t he radio-identified SMG population lies at 
< z >= 2.2 (jChapman et al.ll2005l ). However, the ef- 
fect of the bias toward the true redshift distribution in- 
troduced by the radio selection technique is still under 
debate, as SMGs even beyond z = 4 have been already 
detected by the VLA at 20 cm (e.g..lCapak et al.ll2008l: 
'Coppin et alll2009l: iD addi et al.!i2009a"b!: 'Knudsen et al.l 
2010; Mor rison et al . 2010; Schinnerer et al. 2008). The 
sample of unidentified SMGs could be either spurious 
sources or lie at extreme redshifts but no systematic stud- 
ies on this subsample of SMGs have been conducted so 

far. 

Th e launch of the Herschel observatory (jPilbratt et al.l 
|2010( ) promises a different perspective of SMGs than 
provided by radio observations only and offers a higher 
mapping speed than (sub)mm interferome tric observa- 
tions. Deep PACS (jPoglitsch et al.l l2010f ) imaging at 
100 /xm and 160 /zm on the Great Observatories Origins 
Deep Survey North (GOODS-N) wih sample the FIR 
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emission of these dust-enshrouded high-z objects and 
enable us to study in detail their far-infrared spectral 
energy distributio n (SEP), redshif t distribution, dust 
temperatur es ( Ch anial et alj 120101: lElbaz et all 120101: 
Magnelli et alJ I20l"0l ) and dust masses (jSantini et alJ 



2010( ) . In absence of interferometric observations at mm 



and cm wavelengths, the Herschel-PACS beamsize at 
100 ^m (160 ^m) of 6.7" (11.0") provide a more accurate 
location of the dust emission than the bolometer data 
taken with SCUBA, AzTEC, MAMBO or LABOCA. 
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In this letter we discuss our search for PACS counter- 
parts at 100 /im and 160 /xm, explore the diagnostic 
potential of Herschel-PACS for the counterpart identifi- 
cation and compare it with the widely used identification 
approach using VLA observations. In comparison, the 
Spitzer-MIPS 70 /im and 160 Atm imaging o f the G OODS 
North region performed by iHuvnh et al.l (|2007D in the 
pre-Herschel era detected at relatively high significance 
only two (one) out of 30 SMGs at 70 fim (160 /xm), 
at rather low redshifts (z = 0.5 and z — 1.2). The 
FIR observations presented here will enable us to 
study a significant sample of SMGs in the far-infrared 
wavelength regime. 

2. FAR-INFRARED ASSOCIATION OF SUBMILLIMETER 
GALAXIES 

PACS observations of the GOODS North region at 
100 /im and 160 //m were taken during the Herschel Sci- 
ence Demonstration Phase in autumn 2009 and are part 
of the Guarantee Time extragalactic PACS survey 'PEP: 
The PACS Evolutionary Probe' (PI: D. Lutz). The fi- 
nal images achieve 3cr sensitivities of ~ 3.0 mJy and ~ 
5.7 mJ y at 100 /Ltm and 160 /um respectively. iBerta et"al] 
(|2010r ) describe in their appendix the Herschel-PACS ob- 
servations, data reduction and the blind extraction of 
PACS sources with SNR> 3. Lutz et al. (in prep.) will 
present the fluxes of PACS sources. The PACS images 
are aligned on the GOODS imaging products. Comple- 
mentary to the Herschel PACS observat ions, we use for 
this w ork data from the VLA at 1.4 GHz (iMorrison et al.l 
l2010f) and Spitzer-MIPS at 24 /im (Dickinson et al., in 
prep.). 

In the past years, several groups surveyed the GOODS 
North region using th e bolometer came r as SCUBA, 
AzT EC and MAMBO (iBorvs et al] [200l JGreve et al. 



200g: 'Hug hes et aLlll99a iPerera et al.l 120081: iPope et al. 



2005 : Wang et al.l l2004l)7 These observations discovered 



about 150 SMGs at 850 /im, 1.1 and 1.2 mm. Robustly 
identified VLA and MIPS 24 /tm counterparts are al- 
ready k nown in the literature for SCUBA an d AzTEC 
sour ces (tChapi n et al."200 l iPope era!] 12006'). In addi- 
tion, (Greye^tTal. (2008) presented VLA counterparts for 

11 out of 30 MAMBO sources. 

We search for Herschel-PACS counterparts of SMGs 
that are either detected with a SNR of > 4 or de- 
tected by at least two different surveys. 56 SMGs fulfill 
this criterion which should assure a robust SMC sam- 
ple to work with. Our sample^^ consists of 36 SCUBA, 

12 AzTEC and 8 MAMBO sources and for 15 SMGs 

^* In case that an SMG is detected at subniin and mm bands, 
we give preference the SCUBA detection and include it in our sam- 
ple. For mm-only detected sources we list, if available, the AzTEC 
otherwise the MAMBO detection. 
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Fig. 1. — Top panel: Histogram of distance between the 
submm/mm position from the bolometric map and PACS positions 
of associated FIR sources at 100 ^m (left) and 160 fim (right). 
Black and gray indicate secure and possible PACS counterparts. 
Middle panel: We show the redshift distribution for our SMG sam- 
ple with spectroscopic redshifts (gray) and overplot our PACS de- 
tections at 100 ^m (green) respectively 160 ^m (red). Lower panel: 
We show the SCUBA flux distribution of our SMG sample (gray) 
and overplot our PACS detections at 100 fiva (green) respectively 
160 fj,m (red). 



spectroscopic redshifts (SMG^pec) have been already ob- 
tained. We match our SMG sample with the PACS 
100 /im and 160 /im blind catalogue and search for coun- 
terparts within a radius of 5.5" for MAMBO sources, 
7.5" for SCUBA sources and 9" for AzTEC sources. As 
a sanity check of the blind catalogue, we inspected by 
eye the search region in the PACS imaging. The search 
circles that we have applied correspond to the beam size 
(FWHM) of the different bolometric datasets and should 
guarantee that no reliable associations to SMGs will be 
missed. We also searched the VLA 1.4 GHz map and 
the MIPS 24 /im images for counterparts at the SMG 
positions. 

We calculate the corrected Poissonian probability p 
that an association of SMGs within the search radius 
is a chance coincid ence. This approach (see for details 
iDownes et al.|[l986[) corrects the simple Poissonian prob- 
ability of a detected association for the possibility of 
associations of different nature but similar probability. 
The derived probability of PACS and MIPS association 
is based on raw number counts in GOODS North. We 
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search for VLA counterparts of SMGs down to ia. Ho w- 
ever, the VLA source catalogue (jMorrison et al.l[201Cl( ) is 
only reliable down to 20 fxiy (~ 5cr). Therefore, we assess 
the reliability of VL A counterparts relying on published 
number counts (e.g.. iFtamalont et al.l 120061 ). Similar to 
previous studies, we define following quality criteria for 
assessing the robustness of identified candidate counter- 
parts. We classify association of SMGs with p < 0.05 as 
secure and with 0.05 < p < 0.10 as possible counterparts. 
We uncover PACS secure or possible counterparts for 
22 SMGs, corresponding to a PACS identification rate 
of 39% of our whole SMG sample, see also Table [T] for 
details. Our PACS identificatio n rate of 39% i s lowe r 
than the rate of 54% found by iMagnelU et al.l (|2010D . 
The reason is that'Mag nelli et all ( 2010[ ) focus on already 
radio-identified SCUBA and AzTEC sources with spec- 
troscopic redshifts mainly. The top panel of Fig. [T] dis- 
plays the (sub)millimeter-PACS separation of our associ- 
ation. We did not find evidence of systematic offsets. We 
find PACS counterparts at a mean (sub)millimeter-PACS 

positional offset of A(s„f,mm)mm-PACS100= 5.0" ± 2.1" 

and A(^^^^hmm)mm-PACSwo= 4.9" ± 2.0". This is consis- 
tent with a typical submm/mm position error of about 
3" to 5" and strengthens our choice of the FWHM of 
the bolometric data as search circle for PACS counter- 
parts. The top panel of Fig. [T] shows clearly a dominance 
of secure counterparts. Given this statistical approach, 
associations with offsets much larger than the average 
might not always be correct. 



Table 1 

pacs associations of 56 submillimeter galaxies in 

GOODS-N 

PACS 100 ^m PACS 160 ^tm 

# secure association 

# possible association 
total identification rate 
(sub)mm-PACS offset 



# total association of SMGs 

# only detected at PACS 160 fim 

# only detected at PACS 100 /im 
new CPTs identified by PACS 

# PACS and radio blank fields 



We detect 20 (16) SMGs at 160 (100) ^m and classify 
15 (12) of them as secure and 5 (4) as possible counter- 
parts. The number of PACS counterparts at 160 /im 
is higher than at 100 /im. This is expected as the 
160 /im measurements lie close to the FIR peak. We 
note that based on the corrected Poissonian probability 
p each PACS detection within the bolometer beam (our 
search circle) is classified as associated SMG counter- 
part. Typical PACS fiuxes of these dust-obscured high-z 
sources range between 4.0 mJy to 34.9 mJy at 100 /tm 
and 5.0 mJy to 65.0 mJy at 160 /im. 

We have 43 (28 are secure; 15 are possible) radio- 
identified SMGs in our sample, about 50% of them are 
seen at PACS wavelengths. Vice-versa, only one PACS 
association (secure) is undetected at 1.4 GHz, see de- 
tailed discussion of GN13 (ahas HDF850.4) at end of 
Section [S] One-third of PACS detected SMGs have at 
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least two VLA counterparts. 

3. DIAGNOSTIC POTENTIAL OF PACS OBSERVATIONS OF 
SUBMILLIMETER GALAXIES 

About twice as many SMGs are identified by the VLA 
than by PACS. VLA observations are more sensitive 
for sources at redshifts up to z — 4. None of the 
well-known, spectroscopically identified SMGs at z = 4 
— GN20, GN20.2a/b and GNIO (Daddi et al. 2009a.lj; 
iDannerbauer et al.l2008HWang et al.l2007ll2009( ) — have 
been significantly detected in our PACS imaging, see 
middle panel in Fig. [TJ Another highly promising z > 4 
cand i date. HDF850.1 (ICowie et all l2009t iDunlop et all 
l2004lHughes et al.lll998[ ) is also not seen by PACS. How- 
ever, all these SMGs have radio counterparts. Focusing 
on our spectroscopic subsample of 15 sources, SMGspec, 
we find a trend that the fraction of SMGs detected at 
PACS bands decreases with redshift (see middle panel 
in Fig. [1]), being explained by the fact that PACS fluxes 
drop with increasing redshift. No source beyond z — 2.00 
(GN06) at 100 /tm and z = 2.58 (GN04) at 160 /tm is 
significantly detected by PACS. However, we do not find 
evidence that the PACS detection rate of SMGs corre- 
lates with the SCUBA flux density (Fig. [1]). In addition, 
VLA observations provide subarcsecond accurate posi- 
tions which are essential for the proper identification in 
the optical and near-infrared bands and follow-up spec- 
troscopic observations. PACS cannot deliver positions at 
these accuracies. To summarize, ultra-deep VLA obser- 
vations still remain the best and most effective approach 
to identify SMG counterparts, even up to redshift z — 4. 

Using only radio , 850 /tm and 24 /im flux densities, 
iDaddi et al'T(|2009af) derived accurate {Az/{l + z) - 0.1) 
radio-IR photometric redshifts for SMGs. Motivated 
by this work, we appl y the radio-IR photo metric red- 
shift code presented in iDaddi et al.l ()2009af ) and calcu- 
late photometric redshifts by adding our PACS mea- 
surements to existing SCUBA, MIPS 24 /tm and VLA 
1.4 GHz flux measurements of our SMG sample. Naively 
we would have expected to improve the accuracy of 
the radio-IR photometric redshifts by adding FIR mea- 
surements shortwards to the IR dust peak. In Fig. [2] 
we show the results of o ur attempt using tem plate li- 
braries of local galaxies (jCharv fc Elbajl2001[ ). In the 
left panel w e show photometric redshifts (adapted from 
IDaddi et al . 2009a) based on radio, MIR and submm 
observations for SCUBA counterparts with known spec- 
troscopic redshifts in GOODS North which agree fairly 
well. In the right panel we show photometric redshift es- 
timates adding the PACS measurements to our previous 
measurements in the radio, MIR and submm. Clearly 
the obtained accuracy decreases and these photo-z esti- 
mates are poorer. Omitting certain measurements (e.g., 
MIPS 24/im and/or VLA 1.4 GHz) gives even worse re- 
sults. We obtain similar results by using template li- 
braries from [Dale & Hclou (2003) and th e average SED 
constructed by Michalowski et al] (|201G|) of 73 sp ectro- 
scopically identified SMGs (jChapman et aLll2005l i. This 
result means that most likely the far-infrared proper- 
ties of SMGs are different from templates built to de- 
scribe local galaxies. These observed differences arise 
most probably from the fact that S MGs are colder than 
local ULIRGs of same luminosity (iChanial et all 120101: 
IChapman erall 120051 : lElbaz et al.ll2010t iMagnelh et all 
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I2OIOD . Further investigations are needed to fully exploit 
the Herschel-PACS imaging in order to obtain accurate 
photometric redshifts. 

We already know the spectroscopic redshift for 15 
SMGs in our sample. However, in absence of spectro- 
scopic redshifts for most of our sources, searching for 
correlations between observed properties is essential for 
constraining the nature, redshift and evolution of SMGs. 
We explore if flux ratios involving PACS 100 /xm and 
160 /im measurements can be used as a rough redshift 
indicator. Fig. [3] displays the results of our analysis. 
The PACS color (S'leo /^m/'S'ioo p.m) of SMGs does not 
vary with redshift and is consistent with the predic- 
tion o f models and observed SED s from ICharv fc Elbap 
(120011). iMi chalowski et all (|2010D . lArmus et all (|2007[ ) 
and iPope et al. (2008) . The PACS colors seem to be 
slightly redder than predicted by the various models and 
may indicate a difference between the observed FIR-SED 
of SMGs and templates describing fairly well local in- 
frared galaxies. The composite rest-frame SED for the 
SMGspec (Fig. 131) subsample shows their diversity in the 
far-infrared. Their far-infrared colors are different from 
the template SEDs, most notably for the local luminous 
infrared galaxies. This may explain our finding reported 
above that we are not able to obtain accurate IR-radio 
photometric redshifts by including PACS measurements. 

The PACS fluxes are redshift dependent, whereas the 
(sub)millimeter flux density at z>l is not ( due to the 
negative K-correction; e.g.. iBlain et al.ll2002t ). thus one 
would expect to see a trend between (sub)mm/PACS 
flux ratio for both PACS bands (see also the predic- 
tions based on different templates in Fig. [3]). Indeed, we 
find a fairly clear trend both for Ss5o /xm/'S'ieo tim respec- 
tively 5*850 fim/Sioo fj,m versus spectroscopic redshift for 
the SMGspec subsample-'^^ over the whole redshift range 
spanned by our SMGs from z ~ 1 — 4. We conclude 
that the (sub)mm/PACS flux ratio seems to be a use- 
ful albeit crude redshift indicator and may help to se- 
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Fig. 3. — Left Panels: From top to bottom we plot the flux 

ratios Sieo tim/SloO pm, S'sSO iim/Sieo nm and 5g50 ^m/SloO )im 

versus redshift for SMGs with spectroscopic redshifts. Template 
SEDs of a local LIRG with Lijt Ri 1 x lO" Lq (blue-solid 
line), a local ULIRG with Lju Ri 1 X 10^^ Lq (black-solid line) 
and a HyLIRG with Lju ^ 1 x lO^"^ Lq (red-solid line), all 



from 'Chary & Elbaz? ("2001'), and an average SMG (black-dashed 
line, Michalowski ct al. 2010), are shown for comparison. Average 
SEDs including Spitzer-IRS spectroscopy of local ULIRGs (black- 
dashcd-dottcd line, Armus ct al. 2007) and SMGs (black-dotted 
line, Pope et al. 2008) arc shown as well. The histograms present 
the distribution of SMGs without spectroscopic redshifts (blue- 
flUcd), lower limits of flux ratios are displayed by the pink-empty 
histogram. For the 5i60 f^m/Sioo iim finx ratio we show the upper 
limit as gray-empty histogram. The pink (gray) arrow underlines 
the limits. Right Panels: From top to bottom, we plot the compos- 
ite rest-frame SED — IRAC 3.6 /xm (pink), MIPS 24 /xm (blue), 
PACS 100 ^J.m (green), PACS 160 nm (green), SCUBA 850 /^m 
(gold), VLA 20 cm (red) — of our spectroscopic SMG sample nor- 
malized on SCUBA, VLA and PACS 160 /xm fluxes. Template 
SEDs of ULIRGs and SMGs are shown for comparison. Both data 
points and templates are normalized on an SMG with a SCUBA 
flux of 8 mjy at 2 = 2.2. 




photo-z photo-z 

(Rodlo,850um,24um) (Rodio,850um,PACS 1 60+ 1 00,24um) 

Fig. 2. — Comparison of spectroscopic and photometric redshifts. 
In the left panel we show photometric redshifts based on radio, 
850 ^m and 24 /^m flux densities measurements for SMG coun- 
terpa rts with spectrosco pic redshifts in GOODS North (adopted 
from [Paddi et aP I2009al ) . There is a fairly good agreement. In 
the right panel, we show photometric redshift estimates including 
PACS measurements which are less reliable. 



15 For SCUBA sources from IPope et all II2006I ) we use the de- 
boosted flux. The same we do for AzTEC and MAMBO sources. 
We converted the AzTEC and MAMBO mm flux densities into 
SCUBA flux densities by assuming a flux ratio of Ssso/Si mm of 
2.5 for a ULIRG-like SED at redshift 2-^2-3. 



lect/mark SMGs at very high-redshifts. We explored also 
VLA/PACS and PACS/MIPS 24 ^m flux ratios versus 
redshift but did not observe any correlation, being con- 
sistent with the fact that VLA, MIPS 24 ^m and PACS 
fluxes drop with increasing redshift. Based on this find- 
ing, we investigate the (sub)mm/PACS ratio for the re- 
maining 41 SMGs without any reliable redshift. This 
analysis shows that PACS-detected SMGs tend to have 
lower flux ratios than the PACS-undetected SMG sam- 
ple (Fig. [3]). This trend seems to be more prominent 
for the 5'85o fim/Sieo fxm flux density ratio. Our analysis 
may indicate that indeed as previously discussed PACS- 
undetected SMGs could tend to lie at higher redshifts 
and the (sub)mm/PACS flux ratio could be a crucial tool 
to select the very high-z tail of the SMG population. 

The previous analysis bridges with the discovery that 
12 out of 56 SMGs (21%) in our sample are blank fields 
both at PACS and radio wavelengths. The nature of 
these sources is still not clear. These SMGs could be ei- 
ther spurious sources which is mainly assumed, or at red- 
shifts higher than z > 4. Indeed, how reliable are these 
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Fig. 4. — Prom left to right and top to bottom we show the 50" X 
50" IRAC 3.6/tm, MIPS 24 /^m, PACS lOOAtm and PACS leO^tm 
images of the fie ld of GN13 alias HDF850.4 1 Hugh eTeraLl [l99l : 
IPope et aLl 12009 1, north is up and east to the left. The 7.5" 
search radius circle at the nominal SCUBA position of GN13 (ma- 
genta cross) is drawn. Green (red) circles represent PACS 100/im 
(PACS 160^m) > 3 fT detections and the white square a VLA 
source. IDaddi et al.l |2009a) suggested that the counterpart at 
z = 0.457 in South-West direction ~ 7" away from GN13 is a 
mis-identification. We propose the PACS 100 /im detection on top 
of the nominal SCUBA position (magenta cross) as new counter- 
part at 2 ~ 1 — 2.5. 



(sub)nim sources detections? The sample consists of 
seven submihimeter and five miUimeter selected sources 
(of the latter population, one is seen by AzTEC and four 
by MAMBO). We checked each source for its reliability. 
The mm sources are selected with SNR > 4. Four of 
the se ven SCUBA galaxies ar e listed in both 
(200^ and 'Wang et al.' (2004|V The fact that 
(2006) and Wang et al. (200J) used (nearly) the same 
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SCUBA data, but performed independent data reduction 
and source extraction, is giving more weight to the relia- 
bility of these detections. As a sanity check we inspected 



the SCUBA and MAMBO galaxies in the AzTEC 1 mm 
map. Our analysis suggests that most of these sources are 
reliable. One of them has a SNR ~ 9 at 850 /xm. How- 
ever, we note that only mm interferometric observations 
will reveal unambiguously the reliability of these sources. 
Their non-detection both at radio and PACS wavebands 
in combination with its (sub)mm/PACS flux ratio may 
suggest that these sources could lie at higher redshifts 
than the typical SMC population (e.g., z > 3 — 4). 

We conclude this work by presenting the identification 
of an S MC counterpart havin g a PACS but no VLA de- 
tection. IDaddi et al.l ()2009aD suspected that the source 
at Zg^f^c = 0.457, about 7" away from the SCUBA posi- 
tion, could not be the correct counterpart of GN 13 ahas 
HDF850.4 (iHuehes et al.ll998tlPope"et al.ll2006D . We de- 
tect a faint (4.2cr), secure counterpart at 100 /zm, being 
very close to the bolometer position (~ 0.6"). Its relia- 
bility is fortified by a 200 /uJy strong, secure counterpart 
at 24 /im. Using the flux ratios 5*850 /im/5'ioo /^m = 0.38 
and 5850 ^m/'^ieo /jm > 0.37 as a coarse redshift indica- 
tor suggests a redshift of z ~ 1 — 2.5. 
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